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TECHNICAL PAPER 


TIME-FREQUENCY REPRESENTATION OF A HIGHLY NONSTATIONARY SIGNAL 
VIA THE MODIFIED WIGNER DISTRIBUTION 

I. INTRODUCTION 


This report describes and evaluates a new signal processing technique called the modified 
Wigner distribution (MWD) used for the spectral analysis of highly nonstationary multicomponent 
signals. The regular Wigner distribution (WD) is capable of providing high resolution time-frequency 
estimates of nonstationary signals which are common in many fields including the dynamic response 
of rotating machinery. A traditional tool for such analysis has been the short-time Fourier transform 
(STFT), commonly referred to as “overlap” processing, which is obtained by applying a fixed-length 
moving time window to nonstationary data prior to performing the fast Fourier transform (FFT). 
However, if spectral components within the signal vary considerably in frequency during the time 
window, the STFT often fails in providing enough frequency resolution to identify key time-frequency 
spectral characteristics. 

The WD was originally introduced in 1932 by E. Wigner. 1 It received little attention until 
1980 when Claasen and Mecklenbrauker 2 presented a comprehensive three-part paper describing 
the utility of the WD as a tool for time-frequency analysis of nonstationary signals. 

Major obstacles arise in the direct use of the WD. Most notable of these problems are alias- 
ing and the generation of artifacts or “phantom” spectral peaks in the resultant time-frequency dis- 
tribution. A number of attempts, with varying degrees of success, have been employed to minimize 
these effects. 3-5 This paper presents yet another approach. With the introduction of the “smart win- 
dow,” this approach will hopefully overcome these obstacles. 

The MWD will have important application in turbomachinery diagnostic analysis. It is 
particularly suited to those turbomachinery operations which are highly transient, i.e. during startup 
or ramping conditions. It can also be an important diagnostic tool in failure analysis where the 
dynamic signals are highly nonstationary. 


II. THE WIGNER DISTRIBUTION 


The WD is a powerful tool in determining the time-frequency characteristics of a highly non- 
stationary signal. The technique has been applied successfully in the identification of abnormal 
machine operating states through vibration signal analysis. 3 Other applications of the WD include 
the analysis of time-varying spectra in optics, speech, sonar, and seismic signal processing. 

The WD of a real signal r(t) is given by: 

W(t,f) = J z (r+ t/ 2) z*(f-t/2) ei ln f? dx , 


( 1 ) 



where z(t) is the analytic signal of r(t), z(t) = r(t ) + j i(t), and i(t) is the Hilbert transform (HT) of 
r(t), i(t ) = HT{r(t)}. Let the absence of limits on the integral symbol denote the interval (-««, <*>). In 
equation (1), z*(t) represents the complex conjugate of the analytic signal z(t). Therefore, the WD is 
the Fourier transform (FT) of the product between the original forward signal and corresponding 
backward signal both centered at time t. If the frequency of a subject narrowband signal is either 
monotonically increasing or decreasing within its time window, the product of this forward and 
backward signal will provide frequency cancellation, and, as a result, will estimate a frequency 
corresponding to the signal frequency at the center of the time window. However, the product in 
equation (1), which serves as the kernel of the FT, introduces a frequency summing effect Taking a 
real single component signal with a discrete component at/ 0 as an example, the summing effect 
produces a peak in the WD spectrum at 2/ 0 . For this reason, equation (1) scales the subject time 
signal by two. However, this scaling is only aesthetic and does not relax sampling constraints on the 
time signal. If 2 f 0 is greater than the original Nyquist frequency (determined when the subject time 
signal was originally sampled) aliasing within the WD spectrum will occur. For this reason, an 
analytic signal is used in equation (1) since it has no energy at negative frequencies. Aliasing due to 
the frequency summing effect can now be avoided, and, in addition, frequency difference effects 
(coupling between positive and negative frequency components) are eliminated. Essentially, use of 
the analytic signal in equation (1) returns the sampling rate constraint for a subject time signal back 
to its original Nyquist rate. 

Prior to introducing the MWD, sample graphical illustrations of the WD for cases involving 
simple sinusoids will be presented. In the examples, the WD’s will be processed through convolu- 
tion in the frequency domain since this format conveys the frequency summing and difference effects 
well. 


Before developing case I of the WD, equation (1) will be simplified into a convenient form 
representing convolution in the frequency domain. Taking the original equation and replacing t/ 2 with 
r, gives 


W(t,f) = iz(t+ t) z*(t- t) eJ 2 n f T dt . 


( 2 ) 


Note, with removal of t/2 from the WD equation, / now represents twice the actual frequency of 
spectral components contained in the original time signal. To simplify further evaluation, one should 
just consider one cross section of the WD at t = t 0 

W(t 0 ,f) = \z(t 0 + r) z*(t 0 - t) ei 2 n f x dt . n\ 


Letting 


x(t) = z(t 0 +T ) 
y(t) = z*(t 0 -t) . 

Equation (3) can be viewed as the FT of the product x(t)y(t) 

W(t 0 ,f ) = FT{x{t)y{t)) . 


2 



The FT of the product x(t)y(t ) can also be viewed as a convolution in the frequency domain. Given 
that 

x(t) <=> X(a) 
y(t ) <=> T(«) 

where <=> denotes FT pair. 

Equation (4) can now be written as the convolution 

W(t 0 ,f) = 1 X(a)Y(f-a) da . ( 5 ) 

This form of the WD will be used in the following graphical evaluations. 

WD CASE j : ANALYTIC SIGNAL CONTAINING SINGLE SINUSOID 
For this case, the real signal is a single sinusoid at frequency / 0 

r(t) = cos (2nf 0 t) . 

The imaginary portion of the analytic signal, z(t), is calculated through the HT of the real signal 

i(t) = HT{cos (2 Ttfot)} = sin {2nf 0 t) . 

The analytic signal, z{t), now becomes 

z(t) = cos (2jtf 0 t) +j sin (2 nf 0 t) = e j2n fo‘ . 

The complex conjugate of z(t) is 
z * {t) = e -j2nf 0 t 


and 


z*(-t ) = e j2nfot 

Moreover, the FT of both z(t) and z*(-t) is a delta function shifted by/ 0 . 

FT {ei 2n f ot } = 5(f-f 0 ) 

This sample case is now in proper form with respect to equation (4), with 

x(t) = z{T) , 

y(t) = z*(-t) , (t 0 = 0) 


and 
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X(f) = 8(f-f 0 ) , 

Y(f) = 8(f-f 0 ) . 

With t a = 0, equation (5) can be written as a function of frequency only 

W(f)=lX(cc)Y(f-a) da, (6) 

and the WD for case I can now be developed graphically (fig. 1) through convolution in the frequency 
domain. 

The left-hand side of figure 1 displays the translation of the Y(-a) spectrum during the con- 
volution process defined by equation (6). The right-hand side of the figure shows successive contri- 
butions to the WD spectrum as the frequency, /, varies during the convolution. For this case involv- 
ing the analytic signal of a single sinusoid, a contribution to W(f) is made only when / = 2 f 0 . 
Remembering that / now represents twice the actual frequency of spectral components contained in 
the original time signal, / must be scaled by l /2 to yield a correct frequency value. In this case, the 
true frequency would be 2f 0 /2, or f 0 , which is the expected result. 

WD CASE II: ANALYTIC SIGNAL CONTAINING MULTIPLE SINUSOIDS 

For this case, the real signal is composed of two sinusoids at frequencies of f\ and f 2 

r(t ) = {cos (2jtfit) + cos (2 it fat) } 

This case is identical to WD case I except that the real signal contains two discrete components. 
Using the same method of reduction as in the previous case yields 

X(f) = 5(f-f l )+5(f-f 2 ) , 

Y(f) = 5(f-f l )+5(f-f 2 ) . 

Again, using equation (6), the WD for case II can be developed graphically (fig. 2) through convolu- 
tion in the frequency domain. 

As in figure 1, the left-hand side of figure 2 displays the translation of the Y(-a) spectrum 
during the convolution defined by equation (6), and, the right-hand side shows successive contribu- 
tions to the WD spectrum as the frequency,/, varies. Since the original time signal contains multiple 
discrete components, peaks in the WD spectrum are generated when/=2/i and/= lf 2 . Again, 
scaling by x h yields the correct frequency representations of the spectral components. However, in 
this case, the WD spectrum also exhibits a contribution at/=/i+/ 2 . This contribution is due to cross 
term coupling between the components at/i and/ 2 . This coupling is displayed in figure 2 where the 
peaks b and c of the X(a) spectrum line up with peaks c and b of the Y(\fi+f 2 ] - a) spectrum. This 
cross coupling effect contributes a peak to the WD spectrum at the frequency [/i+/ 2 ]. After scaling, 
this false peak would appear in the WD spectrum at a frequency of |/i+/ 2 ]/2. 

Generation of false peaks due to cross coupling effects severely limits the practical applica- 
tion of the WD since few real world signals are single component. As shown in the previous 
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example, use of the conventional WD in developing time-frequency representations of multi- 
component signals, even with the use of an analytic signal, introduces erroneous spectral com- 
ponents. These “phantom” peaks only confuse the resulting WD spectrum. 


III. THE MODIFIED WIGNER DISTRIBUTION 


In order for the WD to perform as a practical tool in the time-frequency analysis of multi- 
component signals, unwanted erroneous spectral components due to cross coupling among both 
negative, and most notably, positive frequency components in these signals must be eliminated. The 
MWD accomplishes this. Thus, the superior time-frequency resolution of the WD can be attained 
without the generation of “phantom” spectral peaks. 

For a cross section of the time-frequency representation of a real signal r(t), at t = t 0) the 
MWD is introduced and is defined by: 

= i W(a-f/2) X(a) Y(f-cc ) da , (7) 


where 


x(t) = r(t 0 +t) , 
y(t) = r(t 0 -t ) , 


and 

X(f)^x(t) , 

T(/)«yW . 

As equation (7) suggests, the MWD is evaluated in the frequency domain using the FT’s of a real 
time signal centered at t 0 and its respective reversed signal also centered at f 0 .The central trait of 
the MWD which separates the new technique from the traditional WD is its smart frequency 
window function, I1J (/). Use of this window function eliminates the cross coupling of positive fre- 
quency components. This in turn prevents erroneous spectral peaks from entering the MWD time- 
frequency representation during the evaluation of a multicomponent signal. Moreover, use of the 
smart window, IB (/), eliminates cross coupling between positive and negative frequency 
components by preventing their interaction during the evaluation of the MWD. There are some 
limitations on the use of the smart window, and they will be discussed later. The MWD does not 
rely on the use of the analytic signal of the subject time signal; thus, the HT of the original time 
signal is no longer necessary. Finally, since the MWD is evaluated in the frequency domain, aliasing 
due to the frequency summing effect is avoided. 

As with the WD, the evaluation of the MWD will also be illustrated graphically. Where the 
WD was reduced to a more convenient form representing convolution in the frequency domain to 
simplify its graphical presentation, the MWD is actually evaluated through convolution in the fre- 
quency domain. The graphical evaluations of MWD’s which follow include the same two cases which 
were presented for the WD. MWD case I develops the spectrum for a real signal containing a single 
sinusoid while MWD case II evaluates a real signal containing multiple sinusoids. 
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Using equation (7), MWD cases I and II can be developed graphically through convolution in 
the frequency domain. 

MWD CASE I : REAL SIGNAL CONTAINING SINGLE SINUSOID 

The real signal to be evaluated is a single sinusoid at frequency f 0 

r(t) = cos (2 nf 0 t) . 

The required time signals centered at t 0 are 

x(t) = cos {2nf 0 (t 0 +t)) , 
y(t) = cos {2jrf 0 (t 0 -t)} . 

To simplify the evaluation, let t 0 = 0, yielding 

x(t) = cos {2nf 0 (+t)} , 
y(t) = cos {2nf 0 (-t)} . 

Fourier transformation of x{t) and y(t) gives X(f) and Y(f) 

X(f) = 5(f-f 0 )+S(f+f 0 ) , 

Y(f) = S (f+f 0 )+8 (f-f 0 ) . 

X(f) and Y(f) are equivalent with both frequency representations consisting of two delta functions, 
one delta function being shifted along the positive frequency axis by f a and the other shifted along the 
negative frequency axis by f c . Using equation (7), the MWD for case I can be developed graphically 
(fig. 3). 

Figure 3 is essentially the same as figure 1 which displays the development of the WD for a 
signal containing a single discrete component. However, since the MWD operates on real signals, 
the X(a) and Y(oc) spectra also contain negative frequency components. Moreover, figure 3 intro- 
duces the smart window function, Ul (f). The window function is a gate in the frequency domain of 
unity amplitude which translates along the dummy variable (a) axis at one half the rate of the Y(-a) 
spectrum translation during the convolution. The window, as shown in figure 3, has a width in fre- 
quency of 2m. This width parameter is almost inconsequential in the processing of single component 
signals such as in this case, but a proper window width is very critical in the evaluation of multi- 
component signals. As previously stated, the window, ID (/), eliminates cross coupling between 
positive and negative frequency components by preventing their interaction. Moreover, it eliminates 
interaction between components in a multicomponent signal. The window function’s influence on a 
multicomponent signal will be discussed in MWD case II. For a single component signal, unwanted 
cross coupling would have contributed a false spectral peak at zero frequency; however, the gating 
action provided by the window function prohibits this coupling. 

For this case involving a real signal containing a single sinusoid, a contribution to the MWD 
spectrum, M(f), is made only when /= 2 f 0 . Notice, that when the Y(-a) spectrum has translated 
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along the a axis by 2 f 0 , the window function, ID (f) has translated by a frequency of one half of 2 f 0 , 
or ,f 0 . This positioning of the window function permits the desired coupling which only provides a 
contribution to the MWD spectrum at /= 2f 0 . Again, this /must be scaled by */2 to yield a correct 
frequency value. This yields the expected frequency of/ 0 . 

MWD CASE II : REAL SIGNAL CONTAINING MULTIPLE SINUSOIDS 

In this case, the real signal consists of two sinusoids at frequencies of/i and /2 

r(t) = {cos (2k fit) + cos (2nf 2 t)) . 

This case is identical to MWD case I except that the real signal now contains two discrete com- 
ponents. Similar simplification as in the previous case yields the frequency representations of the 
time signals x(t) and y(t) with 

m = 8 (f-f 2 ) + 8 (f-fi) + 8 (f+fi) + S(f+f 2 ) , 

Y(f) = S(f+f 2 ) + 8 (f+f^ + 8 (f-fi) + 8 (f-f 2 ) . 

With these frequency representations, the MWD for case II can be developed graphically (fig. 4) 
using equation (7). 

As in the previous example, figure 4 shows how the positioning of the window function allows 
desired couplings which in turn make correct contributions to M(f). These couplings result in peaks 
at the frequencies 2/i and 2/2 in the final MWD spectrum. However, the key feature of the MWD, 
which separates it from the conventional WD, is its ability to eliminate the cross coupling of positive 
frequency components. This is shown in figure 4 when f =/i+/ 2 . At this value of/, notice that 
the components b and c of the X(a) spectrum line up with peaks c and b, respectively, of the 
Y(\fi+f 2 ]-cx) spectrum. This is a graphic representation of cross coupling between positive 
frequency components. This cross coupling is eliminated by the positioning of the frequency window. 
Notice how the smart window is positioned between the components c and b in the Y([fi+f 2 ]-a) 
spectrum thereby preventing cross coupling. In this situation, the width of the frequency window is 
critical since too wide a window would permit cross coupling. For this reason, the frequency 
resolution of the MWD in differentiating spectral peaks is governed by the window function, ID(/). In 
figure 4, if \f 2 -fi] were less than the window width, 2m, cross coupling between positive frequency 
components would occur and erroneous peaks would be contributed to the MWD spectrum just as in 
WD case II. 


IV. IMPLEMENTATION OF THE MODIFIED WIGNER DISTRIBUTION 

IN DIGITAL FORMAT 


Before developing a digital representation of the MWD, its general form given by equation 
(7) will be simplified. Starting with equation (7) 


(7) 


where 


M(t 0 f) = J Mi(a-fll) X(cc) Y(f-a ) da , 
x(t) = r(t 0 +t) , 
y(t) = r(t 0 -t) , 
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and 


X(f)<=>x(t) , 

Y(f)<*y(t) . 

Again, equation (7) defines a cross section of M(tJ) at t = t 0 . By defining 

P = (fll) - a , 


equation (7) becomes 


= -! U)(-p) X(\f/2] - P) Y(\f/ 2] +P)dp . 

Letting /' = (f/2) yields 

= -SlV(-p) Xif'-P) Y(f'+P) dp . 


Since the window function is symmetric, Ul (-/J) = HI (y5) and equation (9) becomes 

M{t 0 f) = -J U)(p) Xif'-p) Yif'+P) dP . 


( 8 ) 

(9) 

( 10 ) 


Finally, by letting y= -/?, and by invoking symmetry of the frequency window once again, the MWD 
representation becomes 


M(t 0 f) = J W (y) Xif'+y ) Yif'-y) dy . (11) 

Note, that in this form, the MWD represents a frequency which no longer needs scaling. Moreover, 
equation (11) is in a form conducive to digital implementation. The digital form of equation (11) is 

+W 

M(k) = X X(k+i) Y(k-i) , (12) 

i=-m 

where X(k) is the discrete Fourier transform (DFT) of a subject real time series, x(n), centered at 
t 0f and Y(k) is the DFT of yin), the reversal of xin), also centered at t 0 . Mik) represents the 
discrete MWD spectrum of a time signal centered at t 0 , and the summation limits in equation (12) 
constitute the smart window function. Since yin) is the reversal of xin), it can easily be proven that, 

Yik) = X*ik)W N ' k . (13) 

where 

W N = e^ 2n/N) - 


N denotes the number of discrete values in the time history used in calculation of the spectrum with 
n serving as a position index. Finally, by combining equation (13) and (12), the final form of the digi- 
tal representation of the MWD is attained 
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( 14 ) 


+m -(*-/) 

M(k) = 2 X(k+i) X*(k-i ) 

i~-m 

Equation (14) states that the evaluation of the MWD at frequency k is simply the sum of the left- 
hand side and right-hand side of a signal’s FFT spectrum, X(k), with both sides centered at fre- 
quency k, modified by a phase correction term of unity amplitude. Note, that when no window is 
applied (w = 0), the MWD reduces to a special form with an amplitude equal to the power spectral 
density (PSD) but modified by a phase term. 


V. WIGNER AND MODIFIED WIGNER DISTRIBUTION EXAMPLES 


In order to test the capability of the MWD in extreme situations, a simulated sine wave 
varying linearly in frequency at a rate of ~400,000 Hz/s was processed using both the traditional 
STFT method and the MWD. Figure 5 shows the STFT isoplot (logarithmic in amplitude) of a sine 
wave whose frequency first decreases then increases rapidly at this rate within 25 ms. The 
sampling frequency in this case is 10,240 Hz, and the length of the moving time window is 12.5 ms 
(corresponding to 128 discrete data points) with 20 Hanning windows applied. Use of the Hanning 
window makes the effective time window approximately 3 ms. As seen in figure 5, the STFT pro- 
duces a very broad spectral peak since the frequency of the component varies considerably during 
each window. However, figure 6 shows the superior frequency resolution gained by the MWD 
(equation (14)) in processing the same simulated signal. Figure 7 shows the block of raw time data 
used in determining the first spectrum of the isoplots in figures 5 and 6. Along with this raw time 
signal, figure 7 also identifies the actual input simulated frequency at time 0.00625 s, the center of 
the first block of data. This frequency, 2,550 Hz, will serve as the benchmark in comparing the accu- 
racy of the spectra generated by the STFT and MWD. Figure 8 shows the first spectra of figures 5 
and 6 in linear format. The increased frequency resolution attained by the MWD is readily apparent. 
Moreover, the MWD is much more accurate in estimating the actual signal frequency at the center of 
the time window. While the STFT overestimated this frequency by 110 Hz, the MWD estimation 
was within 40 Hz. Again, this increased accuracy can be attributed to the frequency cancellation 
effect inherent to the MWD. Notice in figure 6 that, during the turn in frequency of the single com- 
ponent signal, lobing on the inside turn of the time-frequency representation becomes prominent 
since the frequency cancellation effect inherent to the MWD is limited during this extreme transi- 
tional period in which the signal is neither monotonically increasing or decreasing. 

The isoplot of figure 9 has been developed through the STFT of a simulated sine wave whose 
frequency and amplitude both- varied rapidly within the 400 ms of activity displayed. The length of the 
moving time window in the STFT is 25 ms with a sampling frequency of 10,240 Hz. Within the span 
of each time window, the frequency of the simulated sine wave changes approximately 300 Hz. For 
this reason, STFT of the signal produces a very broad spectral peak. Figure 10 shows the conven- 
tional WD spectrum for the same signal processed using the same parameters as the STFT. The 
increased frequency resolution relative to figure 9 is very evident. Figure 1 1 shows the correspond- 
ing MWD spectrum for the signal, which, in this case, provides just as much frequency resolution as 
the WD. 

As shown in the previous figures, relative to traditional techniques, both the WD and MWD 
provide increased frequency resolution for a single component signal. However, this is not the case 
in the evaluation of multiple components signals. Figure 12 shows the STFT isoplot of a simulated 


9 



signal composed of two sine waves. As before, the frequencies and amplitudes of the discrete com- 
ponents are changing rapidly, and the corresponding spectral peaks are very broad. The WD spec- 
trum for the same signal processed using the same parameters is shown in figure 13. With the multi- 
component signal, the analysis form of the WD introduces an erroneous cross coupling component 
between the two simulated sine waves. This “phantom” peak train confuses the resulting time- 
frequency representation of the multiple component signal. Figure 14 is the corresponding MWD 
spectrum for the multicomponent signal. As can be seen in the figure, a high frequency resolution is 
attained (much improved over that of the STFT method) without unwanted cross coupling terms 
corrupting the time-frequency representation. 

Figure 15 is a PSD isoplot showing the shutdown of an SSME alternate turbopump develop- 
ment (ATD) test rig following a component failure. As seen in the STFT isoplot, following the failure 
which occurs at 615 s into the test, the spectral distribution of the proximity probe signal becomes 
very noisy. In order to study the temporal and spectral characteristics of the signal just prior and 
following the failure, analysis focused on a very short time period around 615 s. Figure 16 shows the 
STFT isoplot for a 300-ms period extending from 615 to 615.3 s. No clear spectral characteristics can 
be identified in this time-frequency representation of the proximity probe measurement. Figure 17 is 
the MWD spectrum for the same period using identical processing parameters. While use of the 
traditional WD would introduce numerous erroneous spectral components in trying to improve upon 
the performance of the STFT, the MWD is successful in providing a much clearer time-frequency 
picture of the time signal without the “phantom” peaks. 

As a final example, figure 18 is an “overlapped” (STFT) isoplot developed during the failure 
investigation of space shuttle main engine u/n 0215. During test 901-666, the engine experienced a 
premature cutoff due to a second-stage turbine-blade failure in the high pressure fuel turbopump 
(HPFTP) u/n 5602R1. The 0-1 -kHz isoplot of figure 18 displays the time-frequency history of the 
fundamental shaft rotational (synchronous) frequency, N, as taken from an external pump-end 
accelerometer. As evidenced by the sudden increase in synchronous amplitude in the figure, the 
pump failure initiated at approximately 3.9 s into the test. The STFT and MWD isoplots for the pump 
speed transducer channel are shown in figure 19. The speed probe is a magnetic-type transducer 
which registers four “blips” with every revolution of the pump shaft, and, for this reason, the 4N 
component dominates both isoplots. Note the increased frequency resolution offered in the MWD 
plot of figure 19, especially during the ramp-down period following the failure. In order to attain more 
insight regarding the failure, attention was focused on the 100 ms timeframe surrounding the failure. 
Figure 20 shows the STFT and MWD isoplots for this extremely short period. Again, notice the 
increased frequency resolution offered by the MWD. Moreover, note the enhancement of the novel 
character of the frequency variation. The frequency separation exhibited in both plots at N, 2 N, 3 N, 
and 5N is currently interpreted as being caused by a phase discontinuity in these respective 
frequency components of the speed signal at the time of the failure. This phase discontinuity of the 
shaft precession can be attributed to the sudden turbine blade loss which instantaneously shifted the 
phase of the driving imbalance force. The frequency branching effect does not appear in the 4 N 
component (which represents the pure shaft rotational motion, or motion other than precession) 
since the speed probe continues to count four “blips” per revolution regardless of the phase 
discontinuity at the time of the failure. It is also important to note that as a result of the branching 
(forking) of the frequencies, an artifact is present in the middle of the fork. This results because the 
difference between the upper and lower frequencies of the branch is within the bandwidth of the 
smart window. 
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VI. CONCLUDING REMARKS 


The modified Wigner distribution has proven very successful in determining time-frequency 
representations of highly nonstationary multicomponent signals. The MWD offers the same time- 
frequency resolution as the traditional Wigner distribution; however, through use of a “smart” 
frequency window, it eliminates the annoying cross coupling artifacts inherent to the WD. This 
ability to eliminate “phantom” spectral components along with its simplistic implementation make 
the modified Wigner distribution a promising real-world signal analysis tool. 
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Figure 1. Graphical evaluation of case I, Wigner distribution. 
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Figure 4. Graphical evaluation of case II, modified Wigner distribution. 
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